function Hv = h(alpha, beta)
k = [cos(alpha)*cos(beta);
    cos(alpha)*sin(beta);
    sin(alpha)];

r{1}=[0 0 0];
r{2}=[1 0 0];
r{3}=[1/2 -0.8660 0];
r{4}=[-1/2 -0.8660 0];
r{5}=[-1 0 0];
r{6}=[-1/2 0.8660 0];
r{7}=[1/2 0.8660 0];


phi = nan(7,1);

for i=1:7
    phi(i)=2*pi*r{i}*k;
end

Hv = exp(1i*phi);
end